Spin transport in a one-dimensional anisotropic Heisenberg model 
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We analytically and numerically study spin transport in a one-dimensional Heisenberg model in 
linear-response regime at infinite temperature. It is shown that as the anisotropy parameter A is 
varied spin transport changes from ballistic for A < 1 to anomalous at the isotropic point A = 1, 
to diffusive for finite A > 1, ending up as a perfect isolator in the Ising limit of infinite A. Using 
perturbation theory for large A a quantitative prediction is made for the dependence of diffusion 
constant on A. 



The one-dimensional spin 1/2 Heisenberg model is one 
of the oldest quantum models [H being also the simplest 
model of interacting quantum particles. Despite being 
exactly solvable by the Bethe ansatz calculating its 
dynamical properties, like transport, is by no means sim- 
ple. Understanding transport in the Heisenberg model is 
important for several reasons. First, it is still not known 
what are the necessary requirements for a system to dis- 
play phenomenological transport laws where a current is 
proportional to the gradient of the driving field. Second 
motivation comes from the condensed matter where one 
would like to understand transport in strongly-correlated 
electron systems. A paradigmatic model is the Hubbard 
model, thought to be related to the problem of high- 
Tc superconductivity, whose low energy excitations can 
be described by the antifcrromegnetic Heisenberg model. 
The Heisenberg model is with a very high accuracy real- 
ized also in the so-called spin-chain materials One- 
dimensional spin chains are receiving increased attention 
also from the cold-atoms community, where they can be 
experimentally realized [4] . Particularly notorious to un- 
derstand is the transport of magnetization. It has been 
actively studied for more than 20 years, yet the subject is 
still hotly debated. In particular, numerical calculations 
sometimes give conflicting results while there are almost 
no rigorous statements, exceptions are a finite Drude 
weight at r = 5] and recently also for T = oo Q, 
both for A < 1, signaling ballistic transport. Various, 
mostly numerical approaches, range from the Mazur's 
inequality 0| , the Bethe ansatz calculation [8] , exact 
diagonalization quantum Monte Carlo .lOj, Lanczos 
method [ll| ^wave-packet evolutio n fl2l . Luttinger liq- 
uid theory 'l3], master equation and correlation 
functions pL5]. Because the model is solvable one would 
be tempted to think that it is ballistic 16]. However, 
recently a solvable diffusive quantum model has been 
found There is also mounting evidence 12l. Il4l. Il5| 
that it is diffusive for A > 1. 

The results presented in the present paper, together 
with other recent works, enable us to give a complete 
picture of spin transport in the linear-response regime in 
the Heisenberg model. For A < 1 the model is ballis- 
tic at infinite temperature Q as well as at any finite or 



zero [H| temperature. At A = 1 our results show anoma- 
lous transport at infinite temperature. It is also plausi- 
ble to expect anomalous behavior at finite temperatures. 
For A > 1 and infinite temperature we show diffusive 
transport. Recent work fl8| shows that as one decreases 
temperature, the diffusion constant increases, possibly 
exponentially fast in 1/T. At temperatures lower than 
the gap transport trivially stops. As one lowers A to- 
wards the isotropic point the diffusion constant diverges 
at any temperature, see also [18]. 

To describe a nonequilibrium situation we couple 
boundary spins of the chain to magnetization reservoirs. 
Time evolution of the density matrix describing the chain 
evolves according to the Lindblad master equation. 



(1) 



where the dissipative linear operator C^^^ describing 
bath is expressed in terms of Lindblad operators L^, 
-C^'Ip) = Ek{iLkP,Ll] + [Lk,pLl]y Reservoirs are 

realized by two Lindblad operators at each end, L\ = 
VT^/]~~ila^ , L2 = VT^I + fiCTi at the left end and 
Lf = Vfy/TT]la+ , Lf = Vt^/T^cj- at the right 
end, cr± = (o-^ ± iay)/2. We fix F = 1. In the diffusive 
regime for A > 1 the diffusion constant does not depend 
neither on the value of F nor on the detailed Lindblad 
operators used to model the bath [l^. Our Lindblad 
operators are such that they always induce a nonequilib- 
rium steady state (NESS) with almost zero average en- 
ergy density, which therefore corresponds to an infinite 
temperature [2^. The Hamiltonian is 
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Small A- In this section we are interested in small 
A, so that the a^a^^i term acts as a perturbation of 
the XX model. Analytical perturbative results will serve 
as a spring-board for the discussion of transport for not 
so small values of A < 1. We split the Lindblad su- 
peroperator into two parts, C = -I- A£(^^\ where 
is the part for A = 0, while £(^^) is the pertur- 
bation, i.e., /;(^^)(p) = i[p,Ej cTjCTj+i]- For A = the 
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NESS solution of this master equation, denoted by po, 
is nondegenerate and ballistic and can be neatly writ- 
ten in a matrix product operator form with matrices of 
fixed dimension 4 2l| . Perturbation series ansatz for the 

NESS for small A is p = po + A ■ pi + A'^ ■ p2 -\ . 

We want to calculate the first two orders in A. Also, in 
all terms we will be interested only in the lowest order 
terms in the driving /x, meaning that we study linear- 
response behavior. Plugging this ansatz into NESS equa- 
tion C{p) = 0, using £*^°)(po) = 0, and then equating 
terms with the same order in A, gets us two matrix equa- 
tions, £(°)(Pi) = -C'-^^Hpo) and £(°Hp2) = ~C^^^\pi). 
A known zeroth-order solution po can be used to get in- 
homogeneous terms in the linear equations for pi, which 
can in turn be used to get p2. The only problem is that 
the number of linear equations is exponentially large in 
the length of the chain n. Nevertheless, several general 
remarks can be made because we know po (it contains 
only crf_„, spin current jk = 2(tT^cr^_^^ - alal^^), and 
their products) and the action of C^^^ and C^^^^: the 
first order term pi does not contain any magnetization 
a J or spin current jk] they appear only in the 2nd or- 
der term p2. Therefore, to the lowest order in perturba- 
tion A the current and the magnetization profile do not 
change! They are the same as for the ballistic XX model 
at A = 0. Because the perturbation of a nondegener- 
ate NESS is nonsingular, fixing n, one can always find 
a sufficiently small A such that the perturbative expan- 
sion will converge and the system is ballistic. However, 
for the transport in the thermodynamic limit the rele- 
vant order of limits is first fixing A and p and only then 
sending n — > (X). This limit is more difficult to treat be- 
cause higher order terms, for instance p2, can grow with 
n faster than the first order term pi , causing the conver- 
gence radius to shrink as n — > cxd. As we will see, this is 
indeed what happens. 

We have found the exact expressions for pi and p2 for 
small n < 7; the coefficients in front of all the terms are 
rational functions with a rather large denominators and 
as such not very transparent. We therefore do not give 
their precise form here but rather focus on their scal- 
ing with n in order to infer the convergence radius of 
the perturbative expansion. Because we are interested 
in the spin transport, we need to know the behavior of 
the correction in p2 that involves and spin current 
jk- It turns out that the correction in the current is the 
same for all sites while the correction to the magneti- 
zation depends on the position. For the spin current, 
the coefiicient in front of the term jfc/2" is , — for 



= 4, 



69235 
513248 



for n — 5, while it is 

56317144998719121983 



45569624481 
243264258368 

for n = 7. If one 



for n 6 and ii7362105703777609136 

looks at the dependence of these coefficients on n one 
notices that it is to a very good approximation linear. 
Fitting gives the dependence c{n) = 0.4286 (n — 2.436), 
with deviations being possibly exponentially small in n. 
Similar corrections, all growing linearly with n are also 
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FIG. 1: Expectation value of the spin current on system size 
n obtained by tDMRG. After sufficiently large ric the current 
converges to a n-independent value, signaling a ballistic spin 
transport. For small A this happens at ric ~ at larger 

the scaling seems to be different. 



found for a'j (the one for is in fact equal to —2 times 
the one for the current, for other spins prefactors are 
larger). To sum up, the expectation value of the spin 
current is to order A^ (and to linear order in p) equal to 
(jk) ~ A* - ^J■A'^ 0.429 (n- 2.436). Because the 2nd order 
correction grows with the system size, the perturbative 
expansion holds only for A < 1 / y/n, and therefore breaks 
down in the thermodynamic limit. Unfortunately, from 
our analytical calculation one therefore can not decide 
about the nature of the spin transport at finite A |22j . 

For larger A we used time-dependent density matrix 
renormalization (tDMRG) simulations |U] with p = 0.02 
to get expectations of magnetization and spin current in 
the NESS. Resuhs are in Fig. [H For A < 1 the cur- 
rent saturates for sufficiently large n, with the saturation 
current monotonically decreasing with A. The system 
is therefore ballistic. For small A < 0.5 the character- 
istic He at which j converges to a constant value scales 
as ^ 1/A^, which is the same as the scaling of the an- 
alytical perturbative result. Even though perturbation 
theory fails as n — t' oo, the scaling A^n — const, appar- 
ently caries over beyond the perturbative result. 

Isotropic At the isotropic point the spin current at 
fixed driving scales as j ~ 1/ -y/n, nicely seen in tDMRG 
data in Fig.[2j The isotropic Heisenberg model at infinite 
temperature therefore display anomalous diffusion, with 
the diffusion constant diverging as D ~ y/n. This is the 
first observation of an anomalous diffusion in a coherent 
(Hamiltonian) quantum system. Furthermore, the mag- 
netization profile along the chain has a nice scaling with 
n and p. As can be seen in Fig. [21 the dimensionless scal- 
ing function looks to be very close to arcsinx, however, 
deviations seen in the figure seem to be larger than the 
finite-size or numerical accuracy effects. 

Large A— Perhaps the most interesting regime is for 
A > 1 where numerical calculations point to a diffu- 
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FIG. 2: Left frame: the scahng of the spin current on the 
system size n for A = 1. The current decays only as 
~ l/y/n (solid line), indicating a superdiffusive transport. 
Right frame: Scaling of the magnetization profile at A = 1 
(two overlapping dashed curves) is very similar to arcsin x 
(red solid curve). For A < 1 the profile is flat (dot-dashed 
curve) . 



sive transport at infinite temperature [14l Il5l | . Because 
analytical treatment seems to be hard, we choose to 
study the case of large A, where the Hamiltonian is 
close to the Ising one. The case of large A can be 
equivalently reformulated with the Hamiltonian H = 

<<y^'^'^+i + '^J^J+i) + with e = l/A [H. 

The NESS for e = is exponentially degenerate. It is 
easy to see that the eigenvector of a dissipative bath part 
is CY^^{l-^lal) = at the left end and £^^*^(1+mO = 
at the last spin. Because of [cr|, 7?e=o] = 0, any operator 
of the form [1 — iA,a\)®x® (I + ^ct^), with x being an ar- 
bitrary combination of \j and cr|, is a zero eigenstate of 
the Lindblad superoperator, i.e., the NESS. There are 
2n-2 independent states of this form. Besides these, 
there are additional NESS states, namely, those for which 
£(e=o)^2;) — Q holds and where x now includes also <y^'^ . 
This increases the degeneracy even further. Because of 
this high degeneracy perturbation theory is more difficult 
than for small A. An important thing to note is that the 
NESS for such Ising chain can support an arbitrary mag- 
netization profile, while the spin current is always zero. 
The Ising spin chain is therefore a perfect insulator. Ex- 
ponentially high degeneracy can now be understood also 
as being due to the isolation of the bulk from the bound- 
aries, so that spins in the bulk "do not know" about the 
reservoirs at the boundaries. High degeneracy is there- 
fore generic and can not be removed by a different choice 
of Lindblad operators. Perturbation e breaks this high 
degeneracy making the NESS nondegenerate. The gap 
between two eigenvalues of the Lindblad superoperator 
with the largest real parts scales as for small perturba- 
tions. This means that if we want to calculate the lowest 
order corrections in the NESS exactly, we have to ex- 
pand it upto order e^. Perturbative expansion can be for 
small jjL written as p = ^^^(l + jipo + /iepi + Me^P2 + •••)• 
The resulting linear equations for unknown po,pi and 
P2 are C^^^po) = 0,C^°Hpi) + J^^'^'^Hpo) = 0, and 
L^^\p2) + C'^^^\pi) = 0, where we have split the su- 
peroperator into an unperturbed part Z!^"^ and the per- 



turbation (p) = I]j CT^crJ+i + crjCTj^_J. Using 

the appropriate ansatz, we have obtained exact solutions 
for small n < 9, however, they are again complicated, 
involving many terms. We only point out features im- 
portant for the spin transport. The first observation is 
that po can contain only terms that are already present 
in the NESS for e = 0. This includes magnetization, 
but not the spin current. Therefore, operators (t| will be 
present in po, while the spin current will be present only 
in the first order term pi (because C'^^^^a'^) will result 
in the current). Spin current is therefore always pro- 
portional to e {— 1/A), while the magnetization scales 
as e''. Writing out the equation involving the coeffi- 
cient c in front of the spin current in pi, one gets 
2c = zi — Z2 — (^i,2'23), where zj. is the coefficient in 
front of CT^ in po and (/i 1,2-23) is the coefficient in front 
of (cicrf "I" ^\'^\)'^% i'^ Pi- Our exact analytical solutions 
for small n's show that the term (/i 1,2 -23) is always equal 
to 2c for n > 4. Therefore, for large n one has an exact 
relation Ac — zi — Z2. This states that if the magneti- 
zation profile is linear on average, then the spin current 
scales as ~ 1/n and the transport is diffusive. Of course, 
showing that zi — Z2 ~ 1/n might be no easier than 
showing this for the current. Exact solutions give the 
expectation value of the current j = /ietr (jfepi)/2" as 
(n- 1) • j7(2pe) = 3, |,2, f|, ^, M, forn = 3, . . . ,8. 
To access the limiting value we have looked at the con- 
vergence of j/eS/fj, ~ {n — l){zi — Z2) with n. If the 
transport is diffusive, this coefficient should converge to 
the diffusion constant D. In the Fig. [3] we plot the val- 
ues of these exact coefficients, together with numerically 
obtained ones for n < 24. The scaling seems to be linear 
in 1 jn enabling us to obtain the limiting value of the co- 
efficient as n — >■ 00. Using this limit we can predict that 
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FIG. 3: Finite-size scaling of the diffusion constant for large 
A. Finite value for n — > 00 indicates diffusion. 

the spin current goes as j x e2p^^ + ©(/ie"^), resulting 
in the diffusion constant 24 1 



Dx2.95/A. (3) 
The dependence of _D on A has been discussed in f25| . 
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FIG. 4: The dependence of the diffusion constant D on A. 
Straight hne is the theoretical prediction (|3]) obtained from 
the perturbative treatment of the Ising model. 

The 2nd order term p2 (as well as the 3rd) does not con- 
tain any corrections to the current. From the analysis of 
the solutions for small n it also does not appear that they 
would grow with n. The convergence radius of the per- 
turbative series is therefore finite and does not decrease 
as n — oo. 

To verify the theoretical prediction for a diffusion con- 
stant ([3]) we have again performed tDMRG simulations 
for a range of A as well as n. For each A a diffusive 
scaling of the current j ^ \/n has been checked and the 
prefactor, being the diffusion constant, determined. The 
results are plotted in FigHl together with Eq.(|31). We can 
see a perfect agreement for large A, with the perturbative 
result ([3]) holding upto quite small A [26| . 

Conclusion- By using perturbation theory in the limit 
of small and large anisotropics A as well as large-scale 
numerical simulations we have shown that for A < 1 the 
Heisenberg model displays ballistic spin transport. At 
the isotropic point transport is anomalous, with the cur- 
rent scaling as ~ Xj^fn. This is the first observation of 
an anomalous transport in a coherent quantum system 
and has strong implications for an unexplained high heat 
conductivity measured in spin-chain materials [27| . For 
A > 1 we show that the transport is diffusive and in- 
versely proportional to A for large anisotropics. Support 
by the Program Pl-0044 and the Grant Jl-2208 of the 
Slovenian Research Agency, the project 57334 by CONA- 
CyT, Mexico, and IN114310 by UNAM is acknowledged. 
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